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Abstract 

If new physics (NP) is present in B — ► tttt decays, it can affect the isospin I = 2 or I = 
channels. In this paper, we discuss various methods for detecting and measuring this NP. The 
techniques have increasing amounts of theoretical hadronic input. If NP is eventually detected in 
B — > 7T7T — there is no evidence for it at present — one will be able to distinguish 1 = 2 and 1 = 0, 
and measure its parameters, using these methods. 
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I. INTRODUCTION 



The decay B d — > 7r + 7r was originally thought to be dominated by a single weak decay 
amplitude, the tree amplitude. As such, it could be used to extract the Cabibbo-Kobayashi- 
MasKawa (CKM) ang.e „ Q. However, it was shown that this decay a. 8 o receives significant 
contributions from the penguin amplitude, which has a different weak phase |2j. This is 
referred to as "penguin pollution." Because of this, a cannot be obtained cleanly from 
measurements of B®(t) — > tt + tt~ . 

However, in 1990 it was shown that it is possible to use isospin to remove the penguin 
pollution and extract a f3j- Briefly, the argument goes as follows. In B — > tttt decays, the 
final 7T7T state can have only isospin / = or / = 2. Thus, there are only two decay paths, 
implying that the amplitudes for the three decays B + — > 7r + 7r°, B® — > 7r + 7r~ and B® — ► 7r 7r 
obey a triangle relation. In terms of the isospin amplitudes Aq (I = 0) and A 2 (1 = 2), the 
amplitudes can be decomposed as [3( 

- V2A(B + -> tt+tt ) = 3A 2 , 

-A(B° d - 7T+7T-) = A 2 + A , 

-V2A(B° d -> ttV) = 2A 2 - A , (1) 

where A and A 2 include both weak and strong phases. The amplitudes for the CP-conjugate 
processes are obtained from the above by changing the sign of the weak phases. From these 
expressions, we see that the triangle relation is \/2A(B + — > 7r + 7r°) = A(B d — > + 
— > 7r°7T ), and similarly for the CP-conjugate decays. 
In the Standard Model (SM), the electroweak-penguin (EWP) contributions to B — > titt 
decays are expected to be very small, so that the amplitude A 2 has a well-defined weak phase: 
A 2 = \A 2 \ exp(«7) exp(i5 2 ), where 7 and S 2 are the weak and strong phases, respectively. In 
the CP-conjugate process, the 1 = 2 contribution is thus simply \A 2 \ exp(— ij) exp(i5 2 ). On 
the other hand, due to the presence of the penguin amplitude, A does not have a well- 
defined weak phase. It can be written A = \A \ exp(i6), where \A \ and 9 are complicated 
functions of both weak and strong phases. The 1 = contribution to the CP-conjugate 
processes is therefore \Aq\ exp(z^), where \Aq\ and 9 differ from \Aq\ and 9. We can now 
count the total number of theoretical parameters. Setting 5 2 = (an overall phase is not 
physical), we find there are six: | ^4.2 1 , \Aq\, \Aq\, 9, 9 and 7. However, there are also six 
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independent experimental measurements |4|: three average branching ratios (-B+o, and 
Boo), two direct CP asymmetries (C+_ and Coo; the direct CP asymmetry C + o is expected 
to vanish in the SM), and one interference CP asymmetry (£+_). (The sub-indices refer to 
the charges of the physical pions in the final state.) We can therefore obtain all theoretical 
parameters, up to discrete ambiguities. In particular, combining the phase of B^-B® mixing 
(/?), we can extract a = n — (3 — 7. 

Obviously, it is also possible to extract the other five theoretical parameters. This was 
shown explicitly by Charles j j . In his analysis, a diagrammatic decomposition of the B — > 
7T7T amplitudes was used: 

- V2A{B + -> 7T+7T ) = T + C , 

-A(B° d ^7r + 7T-) = T + P , 
-V2A(B° d -> ttV) = C-P. (2) 

Here, T and C are (nominally) the color-allowed and color-suppressed tree amplitudes, 
respectively, while P is the penguin contribution. All three contain both weak and strong 
phases. The one complication here is that, while T and C are both proportional to V* b V ud ~ 
exp(ry), the penguin amplitude contains three contributions, corresponding to the flavor of 
the internal quark. These can be written 

PuV: b v ud + p c v; b v cd + p t v; b v td = (p u - p c )v: b v ud + (p t - p c )v; b v td , (3) 

where the (i = u, c, t) contain only strong phases. Above, the unitarity of the CKM 
matrix has been used to eliminate the P c term. The (P u — P c ) term can be absorbed in the 
definitions of T and C, so that 

T = [T tr + (P u - P c )]V* b V ud , C=[C tr -(P u -P c ))V: b V ud , (4) 

where T tr and C tr denote the pure tree-level color-allowed and color-suppressed tree contri- 
butions, respectively. These contain only strong phases; the weak phases have been factored 
out. Thus, the diagrams T and C of Eq. (J2J) actually contain both tree and (P u — P c ) penguin 
contributions, while the penguin diagram P is actually (P t — P c )V^ b V td . 

With these (re) definitions, all three diagrams of Eq. (j2J) have well-defined weak phases - 
the phase of T and C is 7, while that of P is —(3. Thus, in this parametrization, taking into 
account the phase of B^-B® mixing (/?), there are again six parameters: the magnitudes \T\, 
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\C\ and \P\, two relative strong phases, and the weak phase a = tt — [3 — 7. As above, it is 
possible to extract all theoretical parameters. 

Having shown that it is possible to extract a and the remaining theoretical parameters 
from B — ► tttt decays within the SM using an isospin analysis, it is natural to ask how things 
are affected should physics beyond the SM be present. In particular, we want to ascertain 
whether it is possible to detect the new physics (NP) and measure its parameters. If so, we 
want to perform the relevant fits using present data. 

We have noted above that there are two isospin amplitudes in B — > tttt, 1 = and 1 = 2. 
Thus, if NP is present, it can contribute to either of these isospin channels. However, in a 

n 

recent paper |6j, we showed that if no assumptions are made concerning hadronic parameters, 
then the 1 = NP amplitude can never be detected, and only one piece of the 1 = 2 NP 
amplitude can be seen. The main aim of this paper is to investigate what happens if one adds 
hadronic input. In particular, we consider the addition of theoretical assumptions about the 
size of \C/T\ and \P/T\, as well as the strong phases of T and the NP contributions. (The 
former assumptions are on stronger footing than the latter.) 

Note that the decompositions in terms of isospin amplitudes (Aq, A2) and diagrams (T, 
C, P) are equivalent. Throughout this paper we will use a mixed notation, in which the SM 
is described by diagrams, but new physics is described in terms of isospin amplitudes. 

In Section 2, we examine the SM explanation of current B — > tttt data, with no theoretical 
hadronic input. We find a good fit in this case. In Sec. 3, we add theoretical hadronic 
estimates of \C/T\ and \P/T\. In this case, there is a discrepancy in the SM fit, but only 
at the level of about 1.5a. We then consider the addition of new physics. If no theoretical 
input is added (Sec. 4), we again obtain a good fit. We consider the addition of such input 
in Sec. 5. We parametrize the NP using reparametrization invariance. Here we fit for 1 = 2 
and 1 = NP separately. We find that the addition of I = 2 NP gives a fit which is not 
great, but still acceptable. There are not enough observables to measure the 1 = NP 
parameters. In order to make more progress, we change parametrizations in the following 
sections. We assume that the NP strong phases are all small in Sec. 6. We find a good fit 
for 1 = 2 NP, but still cannot fit for 1 = NP. Finally, in Sec. 7, we add a third piece of 
theoretical hadronic input, that the phases of T and the NP are equal. In this case, we are 
able to measure both 1 = 2 and 1 = NP parameters, obtaining good fits in both cases. 
We conclude in Sec. 8. 
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TABLE I: Branching ratios, direct CP asymmetries Cf, and interference CP asymmetries Sf (if 



applicable) for the three B — > tttt decay modes. Data comes from Refs. 



10, 



11 



3; 



averages 







BR[W- 6 } 


Cf 


Sf 


B+ - 




5.5 ±0.6 


0.02 ± 0.07 




B°- 




4.6 ±0.4 


-0.37 ±0.10 


-0.50 ±0.12 


B°- 


-> 7r vr 


1.51 ±0.28 


-0.28 ±0.39 





II. STANDARD MODEL: NO HADRONIC INPUT 



We begin with a review of the standard-model fit to current measurements. In order to 
easily compare with the case of new physics in subsequent sections, our analysis is slightly 
different from what is usually presented. 

As discussed in the introduction, the B — > nix measurements alone allow one to extract 
the CP phase a. However, it is possible to measure the weak phases independently. The 
phase (3 has already been measured very precisely in B®(t) — ► J/ipK s : sin 2(3 = 0.726±0.037 
pj. Note that this value is consistent with SM expectations. If new physics is present in 
B — > ixix decays, it could affect B®-B® mixing. However, the data show that, even if NP is 
present in B®-B® mixing, its effect on (3 is minimal. The measurement of sin 2(3 determines 
2(3 up to a twofold ambiguity. Here we assume that 2(3 ~ 47°, or (3 — 23.5°, in agreement 
with the SM. 

The phase 7 can in principle be measured without penguin pollution through CP violation 
in decays such as B — > DK 8|. Alternatively, 7 can be obtained from a fit to a varietj 
of other measurements, some non-CP-violating. The latest analysis gives 7 = 58.2l| 



anew 

r 0. 



Note that this error includes uncertainties in theoretical quantities. For the purposes of the 
fits, we assume symmetric errors, and take 7 = (58.2 ± 6.0)°. 

It is therefore possible to obtain (3 and 7 (or a = 11 — (3 — 7) independently. Here and 
below, we assume that the SM weak phases are already known, having been obtained from 
other decays. The SM analysis therefore consists of a fit to all B — ► tttt data, along with 
these independent determinations of the CP phases. The latest B — » nix measurements are 
shown in Table |U 
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TABLE II: Results of a fit within the SM. Inputs are the B — ► irir data (Table [IJ and independent 
determinations of the CP phases. All angles are in degrees. 



\T\ 


\C\ 


\P\ 


5 C - 5 T 


5 P — 5 T 


21.8 ±1.1 


18.5 ±1.9 


5.6 ± 1.9 


-66.6 ± 13.6 


-52.9 ± 20.7 



Performed in this way, the fit yields a single solution for the five unknown theoretical 
parameters (Table ITT|). We find xLin/d-o.f. = 0.27/2. We therefore see that the SM can 
acceptably account for the data (leaving aside the somewhat small value of Xmin/ 'd.o.f '.). 



III. STANDARD MODEL: HADRONIC INPUT 



There is one potentially worrisome aspect about the fit described in the previous section: 
naive theoretical estimates of the relative sizes of the magnitudes of the amplitudes yield 
\C tr /T tr \ ~ \{P t V t lV u )/{T tr V* h V ud )\ ~ 0.2 [14]. (T tr and C tr are the tree-level color-allowed 
and color-suppressed amplitudes, respectively.) While the \P/T\ ratio in Table ITT1 is in line 
with these expectations, \C/T\ is quite a bit larger than expected. We therefore see that, 
although the SM can account for the B — > tttt data, it appears to give hadronic parameters 
which are at odds with theoretical expectations. It is therefore worthwhile to re-examine 
the fit, including the theoretical hadronic inputs of \C/T\ and \P/T\. 

Consider first the ratio \C/T\: 



c 




Ctr — {Pit 


-Pa) 


T 




T tr + (Pu 


-Pc) 



(5) 

A numerical value of this ratio depends on estimates of the magnitudes and strong phases 
of Ctr /T tr and (P u — P c )/T tr . The theoretical estimate of \Ct r /T tr \ is ~ 0.2 |l4j . We add a 
theoretical uncertainty and take \Ct r /T tr \ = 0.2 ± 0.1. We also have a theoretical estimate 
of \{P t V t * b V td )/{T tr V* b V ud )\ Q- Including an error, it is taken to be 0.2 ± 0.1. However, 
it is expected that P u and P c are somewhat smaller than P t . We assume no particular 
cancellations, and take \(P U — P c )/T tr \ = 0.1 ±0.1. 

For the strong phases, we proceed as follows. All strong phases are due to rescattering 
from intermediate states Q]. Consider first the SM diagram P c [Eq. (JHJ)]. Its strong phase 
arises principally from the rescattering of the b — > ccd tree diagram, T c . That is, the 
(V - A) x (V - A) T c (b -> ccd) rescatters into the (V - A) x V P c (6 -> uud). The Wilson 



coefficients in the effective Hamiltonian indicate that the size of the rescattered P c amplitude 
is only about 5-10% of that of the T c amplitude. Thus, rescattering costs a factor of about 
10-20. However, the Wilson coefficients also show that the T c diagram is about 10-20 times 
as big as the P c penguin diagram. Thus, the rescattered T c contributes a strong phase of 
0(1) to P c . 

There is also a contribution to the strong phase of P c from rescattering from states 
created by P c itself ( "self- rescattering" ) . However, since rescattering costs a factor of 10-20, 
the strong phase generated from this self-rescattering is much smaller than that generated 
by rescattering from T c . This effect is therefore subdominant, and the strong phase from 
this self-rescattering can be neglected. 

The diagram P u [Eq. (j3J)] also receives a large 0(1) strong phase due to rescattering from 
the b —>■ uud tree diagram, T u . Self-rescattering of P u is negligible. 

We therefore see that the principal source of strong phases is from rescattering from the 
tree diagrams T c and T u . Thus, the strong phases of P u and P c are large and can take 
essentially any value. On the other hand, Pt is essentially real and the strong phase of T tr is 
very small (it is pure self-rescattering). Finally, the strong phase of the diag ram Ctr comes 
from rescattering from T c and is smaller than that of P u /P c , but not negligible. We take it 
to be 0±0.2. 

The ratio \P/T\ is given by 

(P t - P c )V t tV u 



p 




T 





(6) 



[T tr + (p u -p c )}v: b v ud 

The same inputs as above are necessary for the estimate of this ratio. 

In all cases, we assume that the magnitude ratios and strong phases have Gaussian 
and flat distributions, respectively. Allowing all quantities to vary in their allowed ranges 
according to their (assumed) distributions, we find the following estimates: 



0.3 ±0.2 



0.2 ±0.1. (7) 



These are the only two theoretical inputs of ratios of magnitudes of amplitudes, and they 
are thought to be reasonably reliable. In performing the fit in this section we exclude other 
inputs, which have to do with estimates of the sizes of individual diagrams and/or strong 
phases, and which are considered to be less solid. (Note that the above procedure gives 
5 C — 5 T = (0 ± 70)°. This is consistent with the result found in Table ITT1) 



TABLE III: Results of a fit within the SM. Inputs are the B — ► tttt data (Table [J), independent 
determinations of the CP phases, and the theoretical expectations of Eq. Ij7]). All angles are in 
degrees. 



\T\ 


|C| 


|P| 


<5 C — 5 T 


Sp — 5 T 


22.3 ±1.0 


16.7 ±2.0 


5.3 ± 1.5 


-61.6 ± 15.1 


-51.7 ±20.1 



We now redo the SM fit, including the theoretical input of Eq. above. We again find 
a single solution for the five unknown theoretical parameters (Table ITLT]) . The results are 
similar to those found for the fit without theoretical expectations (Table ILT|). However, in 
this case, we have a much worse fit: Xmin/d-o.f. = 6.7/4. 

One explanation is statistical fluctuation. The solution is driven by the unexpectedly large 
branching ratio of B° — > 7t°tt°. As more data is taken, this branching ratio might decrease, 
giving a fit which is more in line with theoretical expectations. (Note that Xmin/ 'd.o.j '. = 
6.7/4 only corresponds to a deviation of about 1.5a. Thus, we must stress that the statistical 
discrepancy in the SM fit is not very strong at present.) A second explanation is that the 
theoretical expectation of \C/T\ is wrong. For example, this can happen if P u is larger than 
expected 3| (if P~c is larger, the size of P will in general be affected). The point is that 
the naive estimates of \Ct r /T tr \ ~ \Pt/T tr \ ~ 0.2 are for the amplitudes uncontaminated by 
penguin contributions, whereas the T, C and P amplitudes which appear in the fit all contain 
P u and/or P c contributions [Eqs. (J3J), (p£J)]. If P u is sizeable, and if the strong phases are 
such that there is destructive interference between the color-allowed tree amplitude T tr and 
P u — P c (for T) and/or constructive interference between the color-suppressed tree amplitude 
C tr and P u — P c (for C), one can account for the large \C/T\ ratio. However, this solution 
is somewhat fine-tuned, and we have no idea why P u should be larger than expected. 

A third explanation, which is the one investigated in this paper, is that new physics is 
manifesting itself in the B — > tttt data. 
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IV. NEW PHYSICS: NO HADRONIC INPUT 



As discussed in the introduction, the B — > titx amplitudes can be written in terms of 
isospin amplitudes or diagrams. In the SM these are related as 

A 2 = ^(T + C)e i \ A 2 = I(T + C)e-^, 
A = \{2T - C)e^ + Pe~ tl3 , A = \{2T - C)e~^ + Pe ip , (8) 

where the diagrams T, C, P include strong phases, but their weak phases have now been 
written explicitly. We therefore see that there are only two areas in which new physics can 
enter: the I = or / = 2 amplitudes. In the most general case, these NP contributions have 
arbitrary phases, so that the isospin amplitudes can now be written 

A 2 = ~(T + C)e* + N 2 , A 2 = l -(T + C)e* + N 2 , 
A = \{2T ~ C)e^ + Pe~^ + N , A = \{2T - C)e^ + Pe i/3 + N , (9) 

where Nq, Nq, N 2 , and N 2 are complex numbers. Although the following discussion is 
general, it is assumed that the NP effects, if present, are roughly the same size as the 
B — > 7T7T penguin amplitudes. There are several reasons for this. First, from the theoretical 
point of view, it is most likely that the NP will affect b dqq loop-level processes. Second, if 
the NP were larger, of tree-level size, it is probable that we would already have seen evidence 
for it, through branching ratios, direct CP asymmetries, etc. Finally, if the NP amplitudes 
are smaller, effects due to the SM EWP amplitudes become important, and the NP will be 
more difficult (if not impossible) to detect. 

Now, it was shown in Ref. |l2j that any complex amplitude can be written in terms of 
two weak phases as follows: 

N = N^e^ + , N = N^e-^ + N^ M e^ M , (10) 

where 

N * A1 = 2i sin (0 A1 - A2 ) ' N * M= 22sin(0 A2 -0 A1 ) ' (U) 
Note that the same complex numbers and iV^ A2 appear in the expressions for iV and 
N. This means that each includes a magnitude and a strong (CP-even) phase, but that the 
weak phases {4>ai an d 4>A2, respectively) have been taken out explicitly. 
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These relations, which are known as reparametrization invariance, can be applied to A^ 
and Nq. The key point is that any two weak phases can be used. We therefore choose to 
express N in terms of 7 and —f3; N 2 is written in terms of 7 and another weak phase. The 
choice of this second weak phase is arbitrary; the only condition is that it be known. For 
convenience, we take it to be zero. Writing N 2 and Nq in this way, one obtains [6( 

A 2 =~(t + c)^ + N 2fi , A 2 = h,t + c)e*> + N 2iQ , 
A = \{2t - c)e^+pe-^ , A = ~(2t - c)e*< +pe i ? , (12) 

where 

t + c = T + C + 3iV 2i7 , 
2t-c = 2T-C + 3iVo, 7 , 

p = P + N ,_ (3 , (13) 



with 



N 2a = i 



N 2 -N 2 
2 sin 7 



N 2 + N 2 .N 2 -N 2 
N 2 , = z- 



2 2tan7 
AT No + N sin/3 .iVo-Ar cos/3 

JV0,7 = ~ ; 7TT- T + I' 



Comparing t 
several things 



2 sin(/3 + 7) 2 sin(/3 + 7) 

N q + Nq sin7 Nq-Nq COS7 

^ " 2 sin (/? + 7) " ' 2 sin(/? + 7 ) ' { ' 

re expressions for the isospin amplitudes in Eqs. (JBJ) and ()12|1. one notices 

First, the expressions for Aq and Aq in these two equations have the same 

form. This implies that one cannot detect 1 = NP without hadronic input. Second, the 

expressions for A 2 and A 2 do not have the same form. Thus, even without hadronic input, 

one can detect NP, but only the iV 2i o piece. Because iV^o is a complex number, there are 

two tests for NP, related to its magnitude and phase (assuming that A^2,o is nonzero). These 

are: a nonzero direct CP asymmetry C + q (if the strong phase of iV^o is different from that 

of t + c) and a difference between the value of 7 extracted from B — > tttt decays and that 

obtained from other, independent measurements. 

Absorbing the NP parameters as above, we can now count the number of parameters. 

There are nine: the four magnitudes of t, c, p and A^o; three relative strong phases, and the 
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TABLE IV: Results of a fit including NP, using the parametrization of Eq. (|12|) . Inputs are the 
B — > tttt data (Table QJ and independent determinations of the CP phases. We have factored 
out the (unphysical) overall phase 5n 20 = a r g-^2,o- The magnitudes are measured in eV and the 
phases in degrees. 



\t\ 


c 


IpI 


1^2,0 1 


St ~ $N 2 


$c ~ <>N 2fi 


dp ~ $N 2 ,o 


6.1 ± 2.7 


9.9 ± 13.7 


12.9 ± 3.2 


9.6 ± 6.5 


81.5 ± 70.5 


-40.5 ± 90.0 


22.3 ± 74.1 


2.8 ± 2.6 


19.8 ± 23.8 


11.4 ± 6.6 


13.2 ± 1.3 


41 ± 108 


-174 ± 9 


-48.6 ± 64.2 


22.8 ± 4.0 


18.2 ± 6.7 


7.3 ± 6.5 


2.7 ± 9.3 


-156 ± 52 


155 ± 32 


157 ± 20 


19.6 ± 3.9 


6.1 ± 22.4 


6.4 ± 1.7 


6.0 ± 8.4 


-19.1 ±43.9 


68.9 ± 174 


-127 ± 35 



two weak phases j3 and 7. Thus, if the values of the weak phases are taken from independent 
determinations, the seven B — > tttt observables are sufficient to perform a fit (this fit is taken 
from Ref. Q.) 

As in the SM fit, we assume that the values of weak CP phases are already known. The 
present B — > mr measurements are detailed in Table |U In fitting to the data, we find four 
solutions, shown in Table ITVl Although the second solution gives a nonzero value of |A^2,o|, 
the other three solutions give values which are consistent with zero. We therefore conclude 
that, if one pays no attention to hadronic quantities, as done in this fit, the present B — » titt 
data exhibit no compelling signs of new physics. This is to be expected since, as shown in 
Table [HI the SM can account for the B — > tttt data if hadronic theoretical expectations are 
ignored. That is, there is no need for NP. 

Note that we obtain xlninl 'd.o.f '. = 0.0049/0, instead of zero as expected. This occurs 
because the current data are slightly inconsistent with the isospin {Ao,^} description. 
That is, for the central values 



2 

3 


^4+0 1 


2 + 




2 _ 


2|A 00 | 


2 


2^1^+01 


IA)I 





cos(5 2 - do) = 31 ^ 7"',, A / UUI = 1.07, (15) 

where \A \ is given by 

W^li-llA^ + lA^+lM'). (10) 



Thus, our fit gives a nonzero value of Xmin/ d.o.f '. 
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V. NEW PHYSICS: HADRONIC INPUT 



In the previous section, we performed a fit without hadronic input. However, as discussed 
earlier, we do have some theoretical hadronic information which can be added to the fit. 
This is done in this section. 



A. Preliminary remarks 

The equations in Eqs. (|T2 ]) -([T1 )) can be inverted to give 

A - A 



p = P + jV ,-/3 
t = T + (N 2 „ + jV , 7 ) 



2i sin(/5 + 7) 



A 2 -A 2 A e-^-A e ip 



2% sin 7 2i sin (/3 + 7) 



c =cH^-N»,,) = -*^ + A f a - A : e :\ 

2«sm7 2«sin(p + 7j 

Nv = (17) 
2« sin 7 

Now, in the preceding sections we introduced three equivalent sets of variables. The set 
closest to experiment is B +0 , B + _, B 00 , C +0 , C + _, C 00 , S' + _. From this we can extract the 
isospin set \A \, \A \, \A 2 \, \A 2 \, S 2 — Sq, Sq — Sq, 62 — So, or, equivalently — c.f. Eqs. (|T7j) - 
the diagrammatic set \t\, \p\, \c\, |A^,o|, S t — S^ 20 , 5 C — Sj^ 20 , and 5 P — Sn 20 . These three sets 
are completely determined from experiment. 

Eqs. f|17|) show us that if theoretical hadronic input is added involving any combination 
of T, C, and/or P, we will be able to probe new physics in iV 2i7 and/or the 1=0 new physics 
amplitudes. In particular, one sees that assumptions about P, T, and C will allow us to 
probe No t -p, N 2yl + Nq }1 , and 2N 2il — No tl , respectively. For example, given the experimental 
fit of \p\ and a theoretical assumption about \P\, N _p is constrained by 

\p\* = |P2| + + 2 |p| CQs6p (lg) 

where 9 P = aig(N _pP*). Thus, \p\ 7^ \P\ is a sign of new physics in N _p. Indeed, 

||p|-|P||<|JV ,-^| < \p\ + \P\, 
\\t\-\T\\ < |iV 2i7 + iVo i7 | < \t\ + \T\, 

12 



\C\ \ < \2N- 



2,7 



0,7 1 



< \c\ + \C\. 



(19) 



We can also identify the combinations that probe N 2n and iVo i7 independently 

\\t + c\ - \T + C\\ < 3|iV 2)7 | < \t + c\ + \T + C\, 
\\2t-c\ - \2T-C\\ < 3|JV 0)7 | < \2t-c\ + \2T-C\. 



(20) 



Recall that one cannot identify new physics in / = without hadronic assumptions. But, 
as we have just shown, it takes only one single theoretical assumption (about \P\) in order 
to be able to probe new physics in / = 0. 

Unfortunately, estimates of the sizes of individual diagrams are thought to be less reliable 
than the two theoretical inputs of ratios of magnitudes of amplitudes in Eq. (J7J). Using 
Eqs. f|T7|) we find 

N ,-f) , N 2)7 + N , 7 p 



P 

T 
C 

T 



P 
t 

c 

7 



T 



T 



2N 2 „-N 0a K- 



T 



t' 

2, 7 + ^"0,7 C 



T 



(21) 



leading to 



N< 



0,-/3 N 2a + A^p, 7 p 
T T t 

2 N 2 ^ - N 0tl jV 2 , 7 + iVo, 7 c 
T + T t 



> 


p 




p 




T 




t 




C 




c 


> 










T 




t 



(22) 



These equations show which combinations of new physics parameters are probed when the 
experimental observable \p/t\ (\c/t\) differs from the SM theoretical prediction for \P/T\ 
(\C/T\, respectively). Sadly, these combinations are rather complicated and, in particular, 
we cannot disentangle N 2n from the 1 = new physics contributions using exclusively 
theoretical predictions for the SM \P/T\ and \C/T\. 



B. Fitting for the new-physics parameters 

We will now consider several scenarios for the NP, attempting to fit for all free parameters 
available within each scenario. If we assume that both 1 = and 1 = 2 NP are present, 
then there are 13 unknown parameters: the seven magnitudes of T, C, P, iVo )7 , No t -p, N 2tl , 
N 2 ,o and the six relative phases (recall that we assume that the SM CP phases are known). 
Given the seven observables, a complete fit would require six theoretical hadronic inputs, 
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which is unmanageable. For this reason, we assume that one type of NP dominates over the 
other, and consider 1 = 2 and 1 = NP individually. 



1. New physics exclusively in I = 2 

Consider only 1 = 2 NP. There are nine unknown parameters: five magnitudes of T, 
C, P, N 2n and A^ 2 , and the four relative phases. With seven B — > -kit observables, two 
additional theoretical hadronic constraints are necessary in order to measure all the 1 = 2 
NP parameters, which we take to be those in Eq. (JJJ). 

Before performing this fit, we can make some general observations. Recall that N 2t o is 
determined from experiment. In the absence of / = new physics, Eqs. (fT3|) become 

t + c = T + C + 3iV 2i7 , 
2t - c = 2T - C , 

P = P , (23) 

which, combined with our estimates for r c and r P , allow us to find explicit analytical formulas 
for T, C, P, and N 2>T We do not include them, since they are not very instructive, but they 
do provide a few interesting bounds, 

T 

1 - 





> 


T 






> 


T 





\p/t\ 

r c - \c/t\ 



2 - c/t 

Moreover, dividing the last two of Eqs. (|2*Hj) . we find 

'2t - c\ 



(24) 



i(8 c -S T ) 



V 



r P e 



i(5 P -5 T ) 



from which 





2t-c 






2 - 




r P 


< r c 




P 







2t-c 



P 



(25) 



(26) 



Given the experimental determination of t, c, and p, this constrains the choices of (r c ,r P ) 
which are consistent with the hypothesis that all new physics appears exclusively in I = 2. 

We now redo the NP fit of Table HVI with the theoretical hadronic input. Table IVI shows 
the result of this fit. We obtain Xmin/d-o.f. = 0.27/0. Note that we expect Xmin = f° r 
zero degrees of freedom (i.e. a solution to the equations). The fact that we don't obtain this 
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TABLE V: Results of a fit assuming only 1 = 2 NP, using the parametrization of Eq. (|12|). Inputs 
are the B — > 7T7t data (Table QJ, independent determinations of the CP phases, and the hadronic 
estimates of Eq. (JJJ). All angles are in degrees. 



\T\ 


|C| 


\P\ 


l^2, 7 | 


\N 2 ,o\ 


5 T - 5n 2 , 




dp — 5n 2 , 






24.0 ±3.8 


7.6 ± 5.6 


5.5 ± 1.7 


10.3 ±11.4 


0.55 ±2.81 


-116 ±72 


-116 ± 195 


163 ± 69 


131 ± 65 





indicates that the equations to be solved are inconsistent, and suggests a poor fit. However, 
at this point we are not certain of how poor this fit is since we cannot test the goodness of 
fit with zero degrees of freedom. 

Fortunately, there is an alternative: consider again the four solutions of Table IIVI which 
correspond to the inclusion of NP without theoretical hadronic input. With hadronic input, 
Eq. (J2T)j) should be satisfied. However, the four solutions in Table ITVl lead to 

= 1.50 ± 0.95, 2.17 ± 3.30, 5.00 ± 6.55, 6.20 ± 1.41 , (27) 

respectively. Using Eq. ([26)1 . we see that none of the central values is consistent with 
(r c , r P ) = (0.3, 0.2). If we relied exclusively on the central values, this would show explicitly 
that none of the solutions in Table ITVl actually provides a good fit, assuming our theoretical 
input is correct. However, the errors are still quite large, so that no such conclusion can be 
drawn. 

There is another way to test the goodness-of-fit: we calculate 

(28) 

Using the values of t, c and p given in Table IIV| we find 

cos(5 T - 8 C ) = {3.33 ± 2.13, 3.25 ± 2.13, 2.57 ± 2.84, 2.10 ± 1.93} . (29) 

All central values are much larger than 1, suggesting a poor fit. However, when we take 
into account the error, the third value is (just) consistent with cos(5 T — 5 C ) < 1, so that this 
value is acceptable. (Note that the third solution also has the smallest x 2 .) 

We therefore conclude that, with this parametrization of the new physics, in which the 
1 = 2 NP has zero degrees of freedom, this type of NP is still permitted, although the fit is 

15 



2t-c 
P 



cos(5 T — 5 C 



4 + rl 
4r n 



4rv 



2t-c 



P 



not perfect. There may be a small deviation of about la. This shows that it is still possible 
to keep (r c ,r P ) = (0.3,0.2) by adding I = 2 NP. 

Note that it is possible to improve the fit by choosing a particular value of the NP weak 
phase. For instance, we could take only N 2a or N 2t o to be nonzero. In both of these cases, 
the number of degrees of freedom is increased to 2, and we obtain acceptable fits. However, 
we prefer not to make any assumptions about the nature of the NP (i.e. its weak phase), so 
we do not consider this possibility here. 

2. New physics exclusively in I = 

The case of / = NP is more complicated. Here there are the same number of parameters 
and constraints — nine — but the effect of the B — > ttti observables is different. First, if 
there is no / = 2 NP, the direct CP asymmetry in charged B decays is automatically zero, 
i.e. C + q = 0. This means that the measurements of BR(B + — > 7r + 7r°) and BR(B~ — > 7r~7r°) 
do not provide constraints on the NP. Instead, they give two different measurements of the 
same quantity, \T + C\. In addition, in the presence of only 1 = NP, the situation is 
analogous to that of the SM: the 1 = 2 amplitude has a well-defined weak phase, while 
the 1 = amplitude does not. In the SM, the 6 B — > tttt measurements are sufficient to 
obtain the various parameters, including the CP phase a. That is, some combination of 
B — > 7T7T measurements is equal to a. However, in our analysis, a is assumed to be known 
independently. Therefore the combination of B — > ixix measurements which is equal to a 
does not provide additional information. 

Thus, for the case in which NP appears exclusively in / = 0, we have only five independent 
B — > ixix constraints, so that we need four additional theoretical inputs to measure all the 
1 = NP parameters. Since we have only the two theoretical inputs of Eq. (J7ji. we conclude 
that it is not possible to measure the 1 = NP parameters. 

We note that one can fit for the 1 = NP parameters if one chooses only one of the two 
terms iVo, 7 or No t -p, in which case there are an equal number of parameters and constraints. 
However, as was the case with 1 = 2 NP, we prefer not to fix the form of the NP at this 
time. 

Notice that the quantities t, c and p can still be determined in this case. The key point 
is that these are observables, derivable from the B — > tttt measurements [Eqs. f)12)l - (jl4)l ] . In 



16 



fact, since we are setting iV^o = here, the values of t, c and p can be determined from the 
SM fit in Table ITT1 with the replacements T — > t, C —>■ c, and P — > p. What we cannot do is 
to disentangle all the parameters T, C, P, iV 0j7 , and A^ ,-/3 from our knowledge of t, c and 
P- 

Fortunately, although in this case it is impossible to fit all parameters, we may still gain 
some information about the new physics. In the absence of / = 2 new physics, Eqs. (|13|) 
become 



t + c 
2t-c 
P 



T + C , 

2T-C + 3iVo, 7 , 



(30) 



while Eqs. (J2TJ) become 
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The last equation leads to 



r c - |c/t| 



1 + c/t 



< 



< 



N 0tl p 
T t' 

c 
t 



r c + \c/t\ 
\l + c/t\ 



(31) 



(32) 



meaning that, if we are lucky, a hadronic input on r c might be enough to detect the presence 
of iVo !7 new physics. Using the central values for r c and for t, c, and 5 C — 5 t from Table HTl 
we find 



0.35 < 



< 0.74 . 



(33) 



This means that, assuming the scenario in which new physics appears exclusively in the 
1 = channel and that r c ~ 0.3, the conclusion is that there might be evidence for new 
physics in iVo, 7 . 

Given the bound on \Nq j7 /T\ in Eq. (|32|). one may check whether the magnitudes r P , 
\p/t\, and \(N 0jJ /T) (p/t)\ can close a triangle. If they cannot, then the first of Eq. (j3*Tj) 
can be used to place a lower bound on \N _p/T\, thus detecting new physics in No -p. We 
conclude that, given appropriate experimental numbers and the two hadronic inputs r c and 



one might be able to get information on both \N Q ^/T\ and \N 0t _p/T\. However, the 



current experimental data, which implies that |iVo )7 | 7^ in this scenario — c.f. Eq. (I33|) 
is consistent with N _p = 0. 
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Summarizing, we see that the theory input (r c ,r P ) ~ (0.3,0.2), along with the current 
experimental results, can be accommodated by NP in the 1 = channel. (Note that the 
errors are still very large, so that new physics is not required.) 

VI. NEW PHYSICS: STRONG PHASES 

In the previous section, we saw that, with theoretical hadronic input, 1 = 2 new physics 
gave a poorish fit, and one could not even fit for all the 1 = NP parameters. In both 
cases, the difficulties could be traced to the large number of theoretical parameters. This 
then begs the question: is it possible to reduce the number of parameters? As we explore in 
this section, the answer is yes, but it requires using a different parametrization to describe 
the new physics. 

In the discussion surrounding the SM fit, we pointed out that strong phases are due 
to rescattering. This has implications for the strong phases of the NP amplitudes. The 
strong phase for a particular NP operator can arise only due to rescattering from itself or 
another NP operator. In either case, we have self-rescattering. However, it is reasonable to 
expect that the strong phases generated by such rescattering are not 0(1), but only about 
5-10%. That is, the NP strong phases are expected to be small [18| . Since we can multiply 
any amplitude by an arbitrary phase, a more accurate statement is that all NP amplitudes 
(both 1 = 2 and 1 = 0) are expected to have the same strong phase. The theoretical error 
incurred by this assumption is expected to be at the level of 5-10%. 

Let us illustrate these features by considering 1 = 0, for example. The 1 = NP 
amplitude is 

No = \N , a \ e<« e <: + \ NoJb \ e <f e <f + (34) 

where "a, b, ..." labels the various possible contributions. Within the above model, the 
strong phases for all contributions to a particular NP amplitude are approximately equal. 
In this case, one can combine all contributions into a single effective NP amplitude, with 
well-defined weak and strong phases: 

iV -|iVo, e/ /|e^ P //e^ P . (35) 

Here, $o,eff * s a wea k phase which changes sign in the CP-conjugate amplitudes; 5 NP is the 
strong phase. The 1 = 2 NP amplitude has a similar expression with the same NP strong 
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phase. 

It is also possible to rewrite the 1 = and 1 = 2 amplitudes using reparametrization 
invariance with the specific choices in Eqs. (|14j). but this is not terribly illuminating, so we 
will stick to Eq. fl35J). 

With this form for the NP amplitude, we write 

- V2A(B + -> tt+tt ) = (T + C)e* 7 + 3 |JV 2 , e //| c*e #? , 

-A(B° d -> tt+tt-) = Te^ + Pe* + \N 2 , eff \ e*%tte«" P + \N 0>eff \ e<°ffe isNP , 
-V2A(B° d -> ttV) = Ce^ - Pe-tf + 2 |iV 2 , e// | e^tfe*"* - |JV ,e//| e 4 <"//e^ P (36) 

The weak phases change sign in the B decay amplitudes. Note that, with this parametriza- 
tion, 1 = 2 and 1 = NP is described by a single amplitude each, instead of two as was 
done in the parametrization of reparametrization invariance. 

As before, we first assume that both 1 = and 1 = 2 NP are present. In this case, there 
are 11 unknown parameters: the five magnitudes \T\, \C\, \P\, |-Wo,e//| and \N2, e ff\, four 
relative strong phases, and the two NP weak phases ^oeff an d ®2,eff- Since there are only 
9 constraints (7 B — » ixix measurements, along with two theoretical hadronic inputs), there 
are more theoretical parameters than there are observables, and we cannot perform a fit. 
We therefore assume that one NP amplitude is dominant. 

3. New physics exclusively in I = 2 

For the case in which one has only 1 = 2 NP, there are 8 theoretical parameters and 
9 measurements/constraints. We can therefore perform a fit, whose results are shown in 
Table IVT1 In this case, we find a single solution, with Xmm/d-o.f- = 0.298/1. Note that this 
fit includes the same difficulties with cos(£ T — 8 C ) detailed in Eq. (J2EJ). However, since there 
is one degree of freedom in this case, we obtain a good fit. 

4- New physics exclusively in I = 

For the case with only 1 = new physics, we still cannot perform a fit. As with 1 = 2 NP, 
there are 8 theoretical parameters and 9 measurements/constraints. However, as discussed 
earlier, two of the B — » titt measurements are redundant and do not probe the NP parame- 
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TABLE VI: Results of a fit assuming only 1 = 2 NP, using the parametrization of Eq. (J36|). Inputs 
are the B — > 7T7t data (Table QJ, independent determinations of the CP phases, and the hadronic 
estimates of Eq. (JJJ). All angles are in degrees. We find x ^ /d.o./. = 0.298/1. 



|T| 


|C| 


|P| 




<5 T - <J" P 


5 C - S NP 


5p - (F p 




24.1 ±3.8 


7.5 ±5.6 


5.3 ± 1.4 


10.3 ± 13.2 


114 ±18 


115 ± 189 


33.6 ±42.8 


59.0 ±5.5 



ters. That is, for the purpose of a fit, there are effectively only 7 measurements/constraints. 
We can therefore not measure the NP parameters in this case. 

VII. NEW PHYSICS: MORE HADRONIC INPUT 

On the previous section, we saw that if we assume that the NP amplitudes are described 
by a single term each, we can fit for the 1 = 2 NP, but not for 1 = 0. In the latter case, we 
need an additional theoretical hadronic constraint. 

This comes from the strong phases. We have argued that the strong phase of T is also 
due to self-rescattering, and should therefore be equal to that of NP. We therefore take as 
the third input 

5 T - 5 NP = (0 ± 11)°, (37) 

where I = 0, 2. With this third input, we can perform a fit for all NP quantities, in particular 
the 1 = NP parameters. 

In this section we perform fits to both 1 = 2 and 1 = NP. In order to obtain the 1 = 
NP parameters, we need all three theoretical inputs. For 1 = 2 NP, we need only one, but 
for simplicity we use all three. 

For 1 = 2 NP, the fit yields one solution for |A^2, e //| (Table IVIT]) . This solution is rather 
intriguing, as it yields a good value of Xmin/ d-o.f '. and a nonzero value of \N 2te ff\ at more 
than 3a. However, we must emphasize that this conclusion must be viewed with a great 
deal of skepticism since this fit has much theoretical input. Indeed, though we did our best 
to estimate the theoretical quantities, we could easily have incorrect central values and/or 
errors. Still, the possibility of NP in B — ► tttt decays is rather interesting. As we noted in the 
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TABLE VII: Results of a fit assuming only 1 = 2 NP, using the parametrization of Eq. Q36JI . Inputs 
are the B — > 7T7t data (Table QJ, independent determinations of the CP phases, and the hadronic 
estimates of Eqs. Q and (37|). We find xL»/d.o.f. = 0.66/2. 



|T| 
<5 T - <J" P 


|C| 
<5 C - 5 NP 


|P| 

5 P - S NP 


W2,eff\ 


24.8 ±4.1 


8.1 ±5.8 


5.9 ±1.3 


7.0 ±2.0 


-3.7 ± 7.5 


10.5 ±23.7 


-131 ± 17 


58.9 ±5.7 



TABLE VIII: Results of a fit assuming only 1 = NP, using the parametrization of Eq. (|36|) . 
Inputs are the B — ► 7T7t data (Table [J), independent determinations of the CP phases, and the 
hadronic estimates of Eqs. and (|57|). We find xL in /d.o.f. = 0.84/2. 



|T| 
<5 T - <F P 


\C\ 

6 C - 5 NP 


|P| 

5 P - 5 NP 


\N ,eff\ 


40.0 ± 18.4 


16.4 ±9.3 


8.0 ±5.4 


20.4 ± 17.8 


4.6 ±9.9 


-118 ±54.5 


-28.9 ±28.9 


-86.6 ±23.1 


40.0 ± 18.4 


16.4 ±9.3 


8.0 ±5.4 


20.4 ± 17.8 


4.6 ±9.9 


-118 ±54.5 


-28.9 ±28.9 


-128 ± 18.7 



discussion of the SM (Sees. 2 and 3), the large value of \C/T\ in the SM fit may be pointing 
to the presence of new physics. Furthermore, the latest measurements of B — ■> 7rK seem to 
indicate the presence of NP in that sector 19]. The two NP solutions may be related. 

In the fit for / = NP, there are two solutions for |iVo, e //| (Table IVIII|) . In this case, 
because the number of degrees of freedom is two, both fits are acceptable. 



VIII. CONCLUSIONS 



If there is new physics (NP) in B — ► tttt decays, it can affect the isospin / = or / = 2 
modes. In this paper, we have discussed several ways of measuring these NP parameters. If 
no hadronic input is used, one can never detect 1 = NP, and only one piece of the 1 = 2 
NP can be observed. Thus, in order to fully detect and measure the NP, it is necessary to 
add theoretical hadronic input. 
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The most obvious input consists of theoretical evaluations of the amplitude ratios \C/T\ 
and \P/T\, where T and C are the color-allowed and color-suppressed tree amplitudes, 
respectively, and P is the penguin contribution. (Note that the C and T amplitudes that 
appear in the B — > tttt also contain penguin contributions.) Given this input, there is 
enough information in principle to measure the 1 = 2 NP parameters, but not those for 
1 = 0. In practice, present data yields a fit for I = 2 NP which is not great, so that, with 
this parametrization of new physics, in which the 1 = 2 NP has zero degrees of freedom, this 
type of NP is disfavored (but not ruled out). For 1 = NP, although one cannot perform 
a fit to measure its parameters, one can compute limits to detect its presence. The current 
data suggest that this type of NP might be present, but the errors are still very large. 

One can then make a well-motivated assumption that the NP strong phases are all about 
the same size. This leads to a new parametrization of the NP with more degrees of freedom. 
In this case, the addition of / = 2 NP yields a good fit (although the presence of NP is not 
favored), but one can still not fit for all the 1 = NP parameters. 

Finally, it is then possible to add a third reasonable assumption, that the strong phases 
of T and the NP are roughly equal. In this case, one can again fit for 1 = 2 NP, obtaining 
a good fit. More importantly, there are now enough constraints to measure the 1 = NP 
parameters, and one obtains a good fit in this case as well. 

We must stress that there is no compelling evidence for NP in the current B — > tttt data 
- the SM can account for it with only a 1.5a discrepancy. However, should statistically- 
significant evidence for NP be found in B — > tttt decays, the various methods discussed 
in this paper will allow us to distinguish between 1 = 2 and 1 = NP, and measure its 
parameters. 
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